Non-Local Product Rules for Percolation 



Saulo D. S. RcislB Andre A. Moreira0 and Jose S. Andradc Jrlll 
Departamento de Fisica, Universidade Federal do Ceard, 60451-970 Fortaleza, Ceard, Brazil 

(Dated; December 5, 2011) 

Despite original claims of a first-order transition in the product rule model proposed by Achlioptas 
et al. [Science 323, 1453 (2009)], recent studies indicate that this percolation model, in fact, displays 
a continuous transition. The distinctive scaling properties of the model at criticality, however, 
strongly suggest that it should belong to a different universality class than ordinary percolation. Here 
we introduce a generalization of the product rule that reveals the effect of non-locality on the critical 
behavior of the percolation process. Precisely, pairs of unoccupied bonds are chosen according to a 
probability that decays as a power-law of their Manhattan distance, and only that bond connecting 
clusters whose product of their sizes is the smallest, becomes occupied. Interestingly, our results for 
two-dimensional lattices at criticality shows that the power-law exponent of the product rule has 
a significant influence on the finite-size scaling exponents for the spanning cluster, the conducting 
backbone, and the cutting bonds of the system. In all three cases, we observe a continuous variation 
from ordinary to (non-local) explosive percolation exponents. 

PACS numbers: 64.60.ah, 64.60.al, 89. 75. Da 



The percolation paradigm represents a formidable ex- 
ample where a simple geometrical construction leads to 
profound concepts in statistical physics, with special em- 
phasis on phase transitions, and real applications in sci- 
ence and technology [H-O- Standard percolation pro- 
cesses arc based on local rules, since they are accom- 
plished through random allocation of sites or bonds, 
therefore disregarding any spatial correlation or global in- 
formation involved in the occupation of other elements on 
the lattice. However, in the case of long-range spatially- 
correlated percolation [3, 0], the probability for a site 
to be occupied depends on the occupancy of other sites. 
Moreover, it has been shown that spatial long-range cor- 
relations in site occupancy can give rise to important 
changes on the structural characteristics of the spanning 
cluster as well as its corresponding conducting backbone 
Q . These changes are strong enough to modify the scal- 
ing exponents of traditional (local) percolation. 

Recently, a new percolation model has been proposed, 
the so-called Product Rule (PR) percolation, in terms of 
a bond occupation process that is essentially non-local 
0. In this model, at each step, two unoccupied bonds 
are randomly chosen and associated with weights given 
by the product of the cluster sizes they would poten- 
tially connect. Only that bond which has the smallest 
weight becomes occupied. By comparison with the tra- 
ditional percolation model the PR model presents 
a more abrupt transition when applied to different net- 
work topologies 0-[l3|- potential applications, the PR 
model has been recently associated to the growth dynam- 
ics of Protein Homology Networks [lH as well as to the 
formation of bundles of single- walled nanotubes [lij . 

Regardless of initial claims of a first order transition 
in the PR model 61, h owever, recent analytical and nu- 
merical works [l7. llSj have demonstrated that the al- 
leged "Explosive Percolation" process actually displays a 
continuous, i.e., a second-order phase transition. This 




FIG. 1: (Color online) Pair of unoccupied bonds e; and Cj 
(dashed lines) randomly selected for the application of the 
product rule, according to the probability P{ri) ~ '"7") where 
Tij is the Manhattan distance between sites i and j (black cir- 
cles), and a is a variable exponent. Following the PR model, 
the bond Ci, merging the two clusters in blue (with 3 sites 
each), becomes occupied. The bond Cj would merge the the 
clusters in red (4 sites) and yellow (6 sites), but remains un- 
occupied. 



apparent drawback of the PR model has been some- 
how overstated, in the sense that the model proposed 
by Achlioptas et al. Q certainly represents an original 
and interesting contribution to the field. For instance, 
much less importance has been given to the non-local at- 
tributes of the PR algorithm. As a consequence of this 
non-locality, the model exhibits a percolation transition 
that, although continuous in nature, seems to belong to 
a different universality class than ordinary percolation 
0,[i3,[l3l- Under this framework, the adjacent edge rule 
(AER) model represents a particular case of the PR 
process that is analytically tractable, since the selection 
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FIG. 2: (Color online) Snapshots of the largest cluster at 
Pc for different values of the exponent a, and a lattice size 
L = 64. The bonds forming the conducting backbone are in 
blue, the cutting bones are in red, and the remaining bonds 
of the largest cluster are presented in green. Although no 
major difference can be observed on the mass Mdus of the 
largest cluster, one can notice that the conducting backbone 
occupies a larger fraction of the largest cluster as a increases, 
leading to a substantial decrease on the number of cutting 
bones Ment- 



is restricted to two adjaeent bonds. When applied to 
random graphs, this model still displays a more abrupt 
transition and different scaling properties than ordinary 
percolation. 

In this Letter we introduce a generalization of the PR 
model in which the range of its non-local features can 
be systematically controlled. This is carried out impos- 
ing that pairs of bonds for selection are randomly chosen 
according to a probability that decays as a power-law 
of their Manhattan distance, namely the distance mea- 
sured as the number of connections separating the sites 
in a regular lattice. This physically plausible assumption 
is inspired on a geographical model for complex networks 
where long-ranged shortcuts are incorporated to regular 
lattices. Such a conceptual construction has been ex- 
tensively used as a way to explain the emergence of opti- 
mal navigation and efficient transport in small- world sys- 
tems [i^l ■ As a consequence of the selection rule adopted 
here, we show that the scaling properties of the system 
becomes dependent on the specific value of the corre- 
sponding power- law exponent. A continuous variation is 
then revealed from the traditional to the PR percolation 
behavior. Moreover, the results of our extensive numeri- 
cal simulations provide strong evidence for the fact that 
the AER model, when applied to regular lattices in two- 



dimensions, falls in the same universality class as ordi- 
nary percolation. 

Our bond percolation process takes place on a square 
lattice of size L. At each step, two sites i and j are ran- 



domly selected with probability P{rij) 



where 



is the Manhattan distance between i and measured as 
the number of connections separating these sites in the 
underlying regular lattice [21|. From each site i and j. 



one bond is then selected among its four adjacent edges, 
namely and e^ , respectively. If at least one of these 
two bonds is already occupied, the entire process of se- 
lection is restarted. If not, following the product rule, 
weights are assigned to each of these bonds, in propor- 
tionality to the product of the size (number of sites) of 
the clusters they would potentially connect. In the case 
a bond connects two sites in the same cluster, the weight 
is equal to the square of the cluster size. The bond asso- 
ciated with the smallest weight becomes occupied, while 
the other stays unoccupied, but can be selected again in 
later steps. 

Our model displays two distinct limiting behaviors, de- 
pending on the exponent a. For a = 0, we recover the 
usual PR, for which the preliminary random selection of 
the bonds ei and 62 constitutes a highly non-local pro- 
cess [1, [13 • In the limit of a — )■ 00, the bonds ei and 
62 are always adjacent, which corresponds to the AER 
process proposed in Ref. [3l , but applied here to regular 
lattices. Although more spatially restricted than the PR 
process, the AER is still non-local, since it requires infor- 
mation about the masses of the joining clusters • As 
we show later, the finite low-dimensionality of the square 
lattice employed here attenuates even further the already 
weaker non-local features of the AER process. 

The percolation process stops when one among all clus- 
ters, namely the spanning cluster, connects the lattice 
from top to bottom [ij. At that point, the fraction p of 
occupied bonds corresponds to the percolation threshold 
Pc- For a = 0, we obtain pc = 0.527 ± 0.001, which is in 
good agreement with previous simulation results of the 
PR on the square lattice [ll|, [l^l- Moreover, we find that 
Pc decreases smoothly and monotonically with a from 
this value to 0.522 ± 0.001 at a = 4 (not shown). Next 
we apply the burning algorithm [2^ to compute the mass 
of the spanning cluster Mdus, the mass of its conducting 
backbone Mf,acfc, and the mass (number) Mcut of cutting 
bonds. The last ones, also called red bonds, if removed, 
would break the spanning cluster in two, therefore de- 
stroying the global connectivity of the system. 

Our results show that, regardless of the value of a, all 
these critical quantities scale with the system size L as 
typical power laws, Mtack ^ L''*'"^'', Mcut ^ L'^<=^* (see 
Fig. [31), and Mdus ^ L'^"'"' (not shown), where dback, 
dcut and dcius are the fractal dimensions of the conduct- 
ing backbone, the cutting bonds, and the largest clus- 
ter, respectively. In Figs.ma)-(c) we show that all these 
exponents exhibit a monotonic variation with a, going 
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FIG. 3: (Color online) (a) Log-log dependence of the mass 
of the conducting backbone M^ack on the system size for dif- 
ferent values of the exponent a. (b) the same as in (a), but 
for the number of cutting bonds Mcut- In both cases and for 
all values of a, the evidence of scaling behavior substantiates 
the calculation of the fractal dimensions dback and dcut as 
the slopes of the corresponding straight lines that are best- 
fitted to the simulation data. All quantities are averaged over 
at least 2500 realizations precisely at the point in which the 
largest cluster appears. 



from a saturation regime of (non-local) explosive per- 
colation at a = to another compatible with ordinary 
bond percolation (BP) at sufficiently large values of a. 
Accordingly, for a = 0, we recover the previously nu- 
merically calculated values of ddus = 1^96 ± 0.01 [ol [lo|. 



dback = 1.52 ± 0.03, and dcut = 1.02 [l4|. In aU three 
cases, by increasing a, a crossover from PR to BP takes 
place in the interval 1 < a < 3. More precisely, ddus de- 
creases in this interval and starts fluctuating around 1.89 
for a > 3 (see Fig. Hk) , in agreement with the classical 
2D value of 91/48 [l|. After increasing in the interval 
1 < a < 3, the exponents dbacfc ^md drut remain practi- 
cally constant for a > 3, around 1.64 [2| (see Fig.|3b) and 
0.75 [l[ (see Fig. respectively. These values are fully 
compatible with previously reported numerical calcula- 
tions for ordinary (local) percolation in 2D. The varia- 
tions of the exponents ddus and dback within 1 < a < 3 
reflect relevant changes in compactness of the spanning 
cluster and its conducting backbone. Although the span- 
ning cluster becomes less compact as a increases {ddus 
decreases) , the mass of the backbone M^ack tends to oc- 
cupy a larger fraction of Mdus , since the dimension d^ack 
increases in the same interval of a values, as shown in 
Fig [31 In these circumstances, a more compact conduct- 
ing backbone implies a smaller number of cutting bones 
(see Figure [3]), therefore explaining the decrease in the 
exponent dcut- 

Next we provide some analytical arguments that in- 
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FIG. 4: (Color online) Dependence on the exponent a of the 
size-scaling exponents for (a) the mass of the spanning cluster 
dciua, (b) the mass of the conducting backbone dtacfc, and (c) 
the number of cutting bonds dcut- In all cases, a crossover 
can be observed in the interval 1 < a < 3 from a regime of 
non-local explosive percolation at a = 0, to a regime that 
is compatible with ordinary bond percolation (BP), at suffi- 
ciently large values of a. The dashed red lines correspond to 
ddu/^ 1.96 [3,[3 and 91/48 ^jn (a), dt^^k = 1-52 M and 
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1.64 0] in (b), and dcut = 1.02 [ill and 0.75 0] in (c; 



dicate how non-local features are introduced in our 
percolation model through the power-law probability 
P{rij). We first consider the average distance be- 
tween all pairs of sites in a empty lattice, (r) = 

12r=i '^-^rf ""/ J2r=i ^r?""", whcrc Nr = 4r is the num- 
ber of sites that are at a Manhattan distance r from a 
given site in the square lattice. Approximating the sum 
by an integral, we obtain that (r) ~ r'^~°'dr. It follows 
that, for a < 2, (r) is limited by the network size, lead- 
ing to (r) ^ L, while, for 2 < a < 3, this distance scales 
as (r) ^ L'^~". For a > 3 and sufficiently large lattice 
sizes, (r) is always finite. As a consequence, the effect 
of non-locality on the scaling properties of the system 
would only play a role for a < 3. In addition, distinct 
non-local behaviors should be expected for the intervals 
a < 2 and 2 < a < 3. These characteristics are consis- 
tent with the results displayed in Fig. [4l The observed 
mismatch between expected and numerically calculated 
crossover values of the exponent a is a consequence of 
finite-size scaling effects as well as the fact that the se- 
quential bond allocation leads to the presence of spatial 
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TABLE I: Estimated values of the percolation threshold pc and the scaling exponent ddus for hyper-cubic lattices of dimension 
D calculated using the jump method [l3.[23|. The presented values correspond to averages over a minimum of 2500 realizations 
of systems with sizes up to L = 4096 (-D = 2), 256 (-D = 3), 64 (D = 4), 28 (D = 5), and 16 (0 = 6). 



correlations in the percolation process. These correla- 
tions make the assumption of an ever empty lattice, as 
adopted to compute (r), no longer strictly valid. 

In order to better confirm the role of non-locality on 
the PR process, we perform additional simulations in the 
two limits of the model at higher dimensions, namely for 
(2 = 0, which corresponds to the original PR process, and 
for the AER process, a ^ oo. In these cases, improved 
performance can be achieved by adopting the so-called 
jump method to analyze the behavior of the order pa- 
rameter Mcius 3 2^ ■ For each realization, we compute 
the average fraction p of occupied bonds at which a jump 
takes place, defined as the maximum change on M^us 
from the occupation of a single bond. This value of p cor- 
responds to the percolation threshold pc- The results for 
Pc and dcius in both limits and different dimensions are 
summarized in Table HI Interestingly, the discrepancy be- 
tween the fractal dimensions calculated for PR and AER 
models increases substantially with lattice dimensional- 
ity. Moreover, our calculations suggest that the resem- 
blance between regular BP and the limiting case a ^ oo 
stands only up to five dimensions, when compared with 
previous results from the literature [l|, Q . 

In summary, wc have proposed a generalization of the 
PR model where the range of non-locality in the perco- 
lation process can be explicitly tuned. Our results show 
that this model displays a rich variety of scaling behav- 
iors, going from ordinary to non-local explosive perco- 
lation. We expect our PR model, since it is based on 
a geographical choice of bond pairs, to provide relevant 
physical insights into the role of non-locality on the crit- 
ical properties of percolation. 
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